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We follo-w the Feynman procedure to obtain a path integral formulation of loop 
quantum cosmology starting from the Hilbert space framework. Quantum geometry 
effects modify the -weight associated with each path so that the effective measure 
' on the space of paths is different from that used in the Wheeler-DeWitt theory. 

■ These differences introduce some conceptual subtleties in arriving at the WKB ap- 

^ . proximation. But the approximation is well defined and provides intuition for the 

differences between loop quantum cosmology and the Wheeler-DeWitt theory from 
a path integral perspective. 
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INTRODUCTION 



^ ■ In loop quantum cosmology (LQC) [1, 2] quantum effects become dominant in the Planck 

CS ! era near the big bang and big-crunch singularities. Because of the underlying quantum 
^ ' geometry of loop quantum gravity (LQG), the theory inherits a repulsive force. This force 
is completely negligible -when the curvature is less than, say, 1% of the Planck scale but 
then grows dramatically, overwhelming the classical gravitational attraction and causing a 
quantum bounce that resolves the classical singularity [3-5]. The effect is generic in the sense 
that it holds, for example, in absence or presence of a cosmological constant, anisotropies, 
and spatial curvature [6-11]. From a path integral viewpoint (see, e.g., [12]), on the other 
^ ■ hand, this stark departure from classical solutions seems rather surprising at first. For, in 
the path integral formulation quantum effects usually become important when the action 
is comparable to the Planck's constant h while the action along classical trajectories that 
originate or end in the singularity is generically very large. Thus there is conceptual tension. 
It is therefore desirable to understand LQC results from a path integral perspective. The 
goal of this paper is to fill this gap at the level of rigor that is common in the discussion of 
path integrals. 

LQC is formulated in the Hamiltonian framework -where one -works with a Hilbert space 
and operators, appropriately taking into account the subtleties that arise because of the 
presence of a Hamiltonian constraint and the absence of an external time parameter. In the 
traditional path integral approach, on the other hand, one generally defines the quantum 
theory by integrating over the space of classical geometries, each path being weighted by 
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e**^ where S is the Einstein Hilbert action, supplemented by suitable boundary terms [12]. 
These procedures then leads to the conceptual tension discussed above. To resolve it, we 
will return to the original derivation of path integrals [13] , where Feynman began with the 
expressions of transition amplitudes in the Hamiltonian theory and reformulated them as 
an integral over all kinematically allowed paths. But now there is a twist: We are dealing 
with a constrained system without external time. In this case the analog of the transition 
amplitude is an extraction amplitude — a Green's function which extracts physical quantum 
states from kinematical ones and also provides the physical inner product between them 
(see, e.g., [14, 15, 17]). If the theory can be deparameterized, it inherits a relational time 
variable and then the extraction amplitude can be re-interpreted as a transition amplitude 
with respect to that time [15]. Irrespective of whether this is possible, the basic object 
that encodes quantum dynamics is the extraction amplitude and we are led to start with 
its expression in LQC and cast it as a path integral. The final result will be quite similar 
to the traditional path integrals based on geometrodynamics [12]. However, there are also 
key differences. As one would expect, these emerge from the fact that since LQG is based 
on a quantum theory of geometry, now one cannot simply start with paths represented 
by smooth, classical metrics. It is precisely these differences that lead to large deviations 
from the classical behavior in the Planck regime near the big bang. Thus, it is possible to 
understand the singularity resolution of LQC also in the path integral language provided the 
choice of paths and the weight assigned to them — i.e., the effective measure — appropriately 
incorporate effects of quantum geometry. Our rewriting of the canonical theory in the path 
integral framework and the associated WKB approximation is also of interest in its own 
right as it is likely to provide new insights and tools in LQC. 

In simplest models one can construct two different path integrals from the Hilbert space 
theory. On the one hand, one can choose to express the transition amplitude as an integral 
over paths in configuration space. On the other hand, one can consider an integral over phase 
space paths. In LQC the geometric operator corresponding to the configuration variable has 
a discrete spectrum. As a consequence, if one descends from the Hamiltonian theory, the 
construction of a configuration space path integral leads to a 'sum over discrete histories', as 
in [14-16]. This sum resembles the vertex expansion of spinfoams [18] rather than a standard 
path integral which features continuous paths with weights given by the exponential of a 
phase {i times an action). In this paper, by contrast, we will use phase space paths and 
arrive at a standard form of the path integral. One can then apply usual techniques as 
for ordinary path integrals. In particular, we will be able to carry out the saddle point 
evaluation of the path integral and use WKB methods. 

The paper is organized as follows. We begin in Sec.II by giving a brief summary of the 
LQC model we wish to use. In Sec. HI we derive the phase space path integral for this model. 
Its saddle point approximation is then studied in Sec. IV. We conclude with a summary and 
discussion in Sec.V. In familiar quantum mechanics the WKB approximation provides the 
leading term in an h expansion of the transition amplitude. In Appendix A we recast that 
h expansion to extraction amplitudes of constrained systems. This expansion provides the 
point of departure for the WKB approximation in Sec. IV. Finally, AppendixB recalls some 
technical results from[15] that are used in Sec. IV. 
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II. SOLVABLE LQC 

We will focus on the LQC model that has been analyzed in the most detail [3-5, 21]: 
the k=0, A=0 Priedmann model with a massless scalar field as a source, which has the 
advantage of being exactly solvable [21]. However, it would not be difficult to extend this 
analysis to allow for a non-zero cosmological constant [6, 7], anisotropies [8, 9] or to the 
spatially compact k=l case [11]. 

In the FRW models, one begins by fixing a (spatial) manifold S, topologically M^, Carte- 
sian coordinates thereon, and a fiducial metric given by q°jjdx"'dx^ — dxl +dxl + dxl. 
The physical 3-metric Qab is then determined by a scale factor a; Qab = a^Qab- f^*^'^' ^^^^ 
Hamiltonian analysis one fixes a cubical fiducial cell V whose volume with respect to g^j, is 
Vo so that its physical volume is ^ = a^Vo- The classical gravitational phase space will be 
described by the canonical pair 

where Pa is the conjugate momentum of a; their Poisson bracket is then given by {u, b} = —2. 
The matter phase space is described by the canonical pair (f) and p, with = 1. Note 

two differences from the previous LQC literature [14, 15, 21]: there is no h'j in the definitions 
of 1/ and b and the momentum conjugate to is denoted here by p rather than p(^). The 
first change removes unnecessary complications in the h expansions one often makes while 
working with path integrals while the second just simplifies notation. 

In the quantum theory, the kinematical Hilbert space is a tensor product "Hkm = ^lin^ ^ 
"H^^" of the gravitational and matter Hilbert spaces. Elements ^(z/) of V.^^^ are functions of 

with support on a countable number of points and with finite norm H^'IP := X^i/ 
The matter Hilbert space is the standard one: 'H^j^^ = d(/)). 

Thus, the kinematic quantum states of the model are functions ^'(z/, 0) with finite norm 
II^IP '■— J2u |^(^^, A (generalized) orthonormal basis in 7^ kin is given by \i',(f)) 
with 

= V. <^(0',0). (2.2) 

A notable feature of kinematics of LQC is that the normalization involves a Kronecker delta 
even though z/ is a continuous variable at this stage. Quantum dynamics is encoded in the 
constraint equation, 

- d^{u, 0) = {-f + 0) M/(z/, 0) = 0, (2.3) 
where p = —ihB^ and is the positive, self-adjoint operator acting on H,^^ [22] given by 

e> ^ (VW\ ^\)\ (2.4) 



o 



Here £„ is related to the 'area gap' A = 4\/37r7£pi ^1 = 7^^; where 7 is the Barbero- 
Immirzi parameter of LQC, v acts by multiplication, and e'^^°^ by translation: (e*^"''^) {v) = 
\E'(z/ -|- 2ioh). (Again, there is a small departure from the notation used in the previous 
literature [14, 15, 21] in that was set to A there. This change will facilitate the semi- 
classical considerations.) The explicit form of is the following second order difference 
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operator 



(e*)(i/) = -^[ ^\v{v + UoK) \\v^ 24^1 + 44^) -2z/2^(z/) 
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+ V^kl^^ - 44^)1 k - 24^1 - 44^) • (2.5) 



From (2.5) one can see that the space of solutions to the quantum constraint is naturally 
decomposed into sectors in which the wave functions have support on specific '//-lattices' [4]. 
Furthermore these sectors are preserved by a complete set of physical observables that is of 
direct physical interest. Thus there is superselection and in each superselected sectors the 
configuration variable v assumes discrete values. For definiteness, we will restrict ourselves 
to the lattice v = AnCoh where n is an integer. On this sector v resembles the momentum 
variable of a particle on a circle, whence the conjugate variable h now lies in a bounded 
interval (0, 7r/4)- 

Solutions to the constraint equation, as well as their inner product, can be obtained 
through the group averaging procedure. Given a state I^E^km) in the kinematical space "Hkin, 
a physical state l^phys) (i-e. a solution to the constraint equation) is given by: 

l^phys) = /dae^°^|^kin) (2.6) 

(see, e.g., [15]). Here, we have introduced the h factor in (2.6) for later convenience. (Since 
the dimensions of the constraint are [C] = ML^, a has dimensions of L~^.) A Green's 
function^ for the above transformation is then given by [14, 15] 

A{vf, (f)f; Ui, (t)i) := / da (i//, 0/| eH^^' 0.) ^ (2.7) 
in terms of which (2.6) can be written as 

*phys(i^, 0) = E /d*^' ^(^' '^^ '^O *kin(i^', 0')- (2-8) 
I/' 

In other words, A gives the matrix elements of the 'extractor' that extracts a physical 
state from every (suitably regular) kinematical one. Therefore, it will be referred to as the 
extraction amplitude. 

The inner product between two physical states |$phys) and I^E'phys) is defined as follows. 
Let |$kin) and j^'kin) be kinematical states such that under the extraction map defined by 
Eq. (2.6) they get mapped to the given physical states. The physical inner product is then 
defined by the action of the 'bra' ($phys| on the 'kef ]^kin), or equivalently, 

($phys, *phys) := (^kinl/daei^^l^kin) 

= J] /d0d0'$kin(^/,0)A(//,0;//',0')*ki„(^/',0')- (2.9) 



We will actually restrict ourselves to the 'positive frequency part' as in [14, 15], so that there is an implicit 
0{p) factor in (2.7) where 9 is the unit step function. We do not write it explicitly just to avoid unnecessary 
proliferation of symbols. 
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Thus, in the 'timeless' framework without any deparamctrization, all the information of the 
quantum dynamics is encoded in the extraction amplitude ^(i^/, 0/; i^i, (pi)- We will start by 
finding a path integral expression for this function. 

Remark: Detailed analysis shows that is a viable relational time variable [3-5] and 
one can use it to deparameterize the theory. When this is done, the extraction amplitude 
A{uf,(f)f; can also be regarded as the amplitude for a transition from z/j at 'time' 0^ 

to Vf and 'time' 0/ [15]. Even though the deparameterized theory is closer to more familiar 
path integrals where one also computes transition amplitudes, in the main body of this 
paper we will work in the timeless framework. Although the two are equivalent, the timeless 
framework is technically simpler because it leads to a path integral that directly involves 
matrix elements of O; in the deparameterized framework it would involve matrix elements 
of \/0 which are much more complicated [15]. 

III. PHASE SPACE PATH INTEGRAL 

In ordinary quantum mechanics, path integrals provide an expression for the matrix 

elements of the evolution operator. Feynman [13] first derived it from the canonical theory 
by writing the evolution as a composition of N infinitesimal ones and inserting complete 
basis between these infinitesimal evolution operators. He then arrived at a 'discrete time' 
path integral expression; the continuum path integral was found by taking the limit N ^ oo. 

In the timeless framework there is no 'cvohition' operator — all we have is the constraint 
equation and its solutions — and the extraction amplitude A{vf, (j)f \ z/j, (pi) replaces the tran- 
sition amplitude. Thus, our task is to construct a path integral expression of this object. The 
idea is to mimic the standard Feynman construction [13] but now using the operator e^'*^ 
that appears on the integrand of Eq. (2.7) for the 'evolution' operator, and then performing 
the a integration. 

Let us be more specific. The integrand of (2.7) can be thought of as a matrix element of 
the fictitious evolution operator es"*^. One can regard aC as playing the role of a (purely 

mathematical) Hamiltonian, the evolution time bring unit, i.e. e^"*" = e^*^ with H = aC 
and t — 1. We now decompose this fictitious evolution into N evolutions of length e— 1/N: 

et*^ — Y[n=i^^^^- By inserting complete basis of the form 1 = Z^^, /d0|i^, 0| in 
between each factor we get 

e^"^ |z/i,0i) = ^ /d07v_i...d0i (i/iv,0Ar|et'°'^|i/iv_i,0jv_i)... 

...{iyi,(l)i\ei^"^\uo,(j)o): (3.1) 

where {unAnI = {^fAf \ and |z/o,0o) = \ i'h<Pi)- 

Let us concentrate on the n-th term appearing in (3.1). Notice that since the constraint is 
a sum of two commuting pieces that act separately on 7^^^" and , one has the following 
factorization: 

(i/n+i,0„+i|e^"^|z/„,0„) = (0„+i|e^^"^'|0„)(i/„+i|e-*^"®|i/n). (3.2) 
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The scalar field factor can easily be evaluated by inserting a complete basis in p, 



6„+i|e^'"^^' \(j)n) = /^e*P"(<^"+i-'^")+i^"P". 



(3.3) 



The gravitational factor in (3.2) is less trivial to compute. As usual in the path integral 
construction, we will take A'" 1 (e = 1/N <C 1) and use an expansion in e to compute this 
term: 

(i/„+i|e-i^«®|i/„) = - |eaK+i|0|i/„) + ^(e^), (3.4) 

where the matrix element of © can be obtained from Eq. (2.5) and is given by 



(There arc several equivalent ways of writing this matrix element. Here we chose one that 
is symmetric in z/„ and z/„+i.) As in usual path integral constructions, we now bring-in b, 
the conjugate variable to u. This can be done through the identity 



TT 



(3.6) 



which, when used to rewrite the Kronecker deltas appearing in Eqs. (3.4) and (3.5), leads 
to the following expression for (3.4): 



e 1 



eae 



= —Jo dbn+ie 



-ib„+i(!^„+i-i/„)/2 



TT 



1 — \ea^^ y/VnVn+l 



sin2(46„+i) + C»(e2) 



= —Jo dbn+ie 

TT 



-|6n+i(t'n+i-^'n)/2-iea^|^^i/„fc'„+i 



-■""+■"1+1 „;„2 



sin2(^o6„+i) 



+ 0(e2). 



(3.7) 



By combining (3.2), (3.3) and (3.7), the amplitude (3.1) takes the form 

{lyj, 0/1 e^-^ li/,, 0,) = Yl (^f • • • dhJd<PN-i . . . d0i x 

(^)^/dp^...dpiei^- + 0(e2). 



I/jV-l,...,i^l 



(3.8) 



where 



S 



N 



N-1 
n=0 



0n+l 



2 



e 

SttG 



V^n^n+l '^" "L'^"'^^ Sin^(4^n+l)) (3.9) 



can be heuristically regarded as a 'discrete-time action'. Thus, following Feynman, we have 
obtained an approximate expression of the extraction amplitude as a sum over phase space 
paths, approximation becoming better and better as e shrinks and A^ grows. 

Already, before taking the N ^ oo limit, we can see two important differences from the 
more familiar path integrals. First, whereas in the classical theory the variables u and b can 
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take all real values, our paths in v take only discrete values and the paths in h are bounded; 
the allowed values of v and h are dictated by spectrum of the corresponding operators 
on the superselected sector we began with. Second, we find that the action is not just a 
discretization of the classical action, but includes quantum gravity corrections (namely, the 
so-called 'holonomy type corrections' encoded here in the sin term). What is the origin 
of these differences? After all, we just mimicked the Feynman construction [13] to arrive 
at (3.9) starting from the Hamiltonian theory. Recall however, that Feynman began with 
the standard Schroodinger representation of quantum mechanics. In LQG, because of the 
underlying diffeomorphism invariance the analogous representation is not viable [19]. This 
difference descends to LQC where kinematics is based on a 'polymer' representation that is 
unitarily inequivalent to the Schrodinger representation [20] . In the polymer representation, 
the quantum constraint operator of the Wheeler-DeWitt theory — which would have been 
analogous to the Hamiltonian operator Feynman used — fails to be well-defined. Defining a 
viable constraint operator on the kinematical Hilbert space of LQC requires an appropriate 
incorporation of quantum geometry underlying LQG [5]. And this is directly responsible for 
the two key differences mentioned above. 

The final step in the path integral construction involves taking the hmit N ^ oo. In this 
limit one typically performs the substitutions ej^^o ~^ I ^'^^ (0n+i — <f>n)/^ d(j)/dT, etc, 
to obtain a continuum (in time) action. However, due to the discrete nature of u, it is not 
possible to interpret the (z^n+i — z^rj/e as a derivative. Therefore, we are led to carry out an 
'integrating by parts,' rewriting this term as 



N-l 



e 

n=0 



N-l 



—0 : — ='1^ 



n=0 



+ ^(Mo-6ivi^iv). (3.10) 



Now the formal limit can be carried out at the same level of precision as is common in path 
integrals and one obtains 



A{uf, cPf, u,, cf>,) = Jda J[Du,{t)] [Vbgir)] [Vp{t)] [P0(t)] et^ (3.11) 



where 



Jo 



( p0 + ^ub -afp'- SttGu'^-^^^^ - ^{ufbf - uA) (3.12) 



and the subscript q emphasizes the fact that we have taken the formal limit of a sum that 
included only 'quantum paths' in the geometrical sector. Thus, the integral is over paths 
v{t) taking only the values v — Anigh and the paths 6(t) take values in the range (0, tt/^q). 
Notice that the classically singular paths where b is divergent are excluded, which refiects 
the fact that the quantum dynamics is free of singularities. However, precisely because of 
the unusual restriction on the domain of integration, this path integral is not of the usual 
type. The path integral (3.11) actually resembles the 'sum over histories' of [15], where one 
sums over the same family of i^-paths. Indeed, by performing the bg integral in (3.11) one 
can recover the sum over histories expansion of [15]. 

But it turns out that we can express this path integral also in a more familiar form using 
a clever trick from path integral framework for a particle in a circle [23] . The trick is to use 
the identity 



2tt 

de f{e,m)e''^^ = I dx I de f{e,x)e'''^ (3.13) 
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which holds for any continuous f{0,x) with a 27r period in 9 to convert the discrete sum 
over m to a continuous integral over x. Let us then return to the finite N or 'discrete-time' 
path integral (3.8) and use this identity. Then, we are led to rewrite each sum over i/„ and 
integral over 6„ appearing in (3.8) as 



•''O „ Jo J -oo J -oo 



(3.14) 



(this is done for n — 1,. . . ,N — 1; the integral over remains unchanged). The allowed 
paths now take value over the whole classical phase space: u is no longer restricted to be 
discrete nor is b required to be bounded. As in standard path integral discussions, it is then 
possible to take the formal limit N ^ oo. We obtain^ 

A{uf, <l>f; u„ (t>i) = !daj[Vv{T)] [^(t)] [Vp{t)] [P0(t)] et^ (3.15) 

where 

[\T(p<j>--bu-a(p^-37rGu^^-^^^^]]. (3.16) 
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While this is the same action as before, one now integrates over all trajectories in the clas- 
sical phase space as in usual path integrals. In particular, the trajectories that make the 
action stationary lie in the domain of integration (which, by contrast, were not included 
in the 'quantum paths' of the previous path integral). We can therefore use all the stan- 
dard techniques such as the saddle point approximation. Thus, while they are (formally) 
equivalent, this second form (3.15) of the path integral is much more convenient in practice, 
particularly to address the issues we began with in Sec.I. This is the path integral form of 
the extraction amplitude we were seeking. 

Since we now integrate over all paths in the classical phase space, in particular, the paths 
are allowed to go through the classical singularity. How can then we see the singularity 
resolution in this setting? The answer is that the paths are not weighted by the standard 
FRW action but by a 'polymerized' version of it which still retains the memory of the 
quantum geometry underlying the Hamiltonian theory. As we will see, this action is such that 
a path going through the classical singularity has negligible contribution whereas bouncing 
trajectories give the dominant contribution. 

Remark: There are other systems in which the passage from the Hamiltonian quantum 
theory to a path integral results in an action that has ^-corrections. Perhaps the simplest 
example is that of a non-relativistic particle on a curved Riemannian manifold for which 
the standard Hamiltonian operator is simply H — —{l^/2m)g'^Va^h- Quantum dynamics 
generated by this H can be recast in the path integral form following the Feynman procedure 
[13]. The transition amplitude is then given by [25] 

(g,%',t') = Pfe(r)]ei^ (3.17) 



^ Note that the integral over a is an ordinary one variable integral. It can nevertheless we reinterpreted as 
an integral over all possible values a(r) together with the a gauge fixing condition da/dr = 0. Doing so 
allows one to rewrite the path integral with any other gauge fixing condition; see for instance [24]. 
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with 

S^JdT{fg,,q-q'^ + £^R) (3.18) 

where R is the scalar curvature of the metric Qab- Extrema of this action are not the geodesies 
one obtains in the classical theory but rather particle trajectories in a ^-dependent potential. 
The two sets of trajectories can be qualitative different. 

IV. SADDLE POINT APPROXIMATION 

As we saw in Sec. II, the extraction amplitude encodes the entire content of quantum 
dynamics. However, in practice, it is difficult to work with the series (3.8) or evaluate its limit 
(3.15). In quantum mechanics and quantum field theory the steepest descent approximation 
is a powerful practical tool to calculate leading contributions to the transition amplitude in 
an h expansion. In particular, this approximation provides the much needed intuition on 
when quantum corrections are dynamically important and when they are not. In Appendix 
A we recast this h expansion in a form suitable for the extraction amplitude of constrained 
systems. We will now use those results to obtain the leading term using a saddle point 
approximation. 

In this approximation, the extraction amplitude (3.15) is given by 

A{uf,cPf;ui,cP,) ^ (det^^^lo)"'^' e^^l" . (4.1) 

Here 5*10 is the action evaluated along the trajectory extremizing the action with initial and 
final configuration points fixed. For given initial and final points, there exist in general two 
trajectories joining them, one with positive and the other with negative p values. As pointed 
out in Sec.II, following the logic spelled out in [15] here we restrict ourselves to the 'positive 
frequency' branch, and so only the p > trajectory gets picked. We will evaluate the phase 
factor in Sec. IV A. The prefactor (det^^S'lo) represents a formal infinite dimensional 
determinant which we will evaluate in Sec. IV B. In Sec. IV C we compare the resulting 
approximate amplitude to the exact one, computed numerically. 

Before proceeding with these calculations, we would like to point out a conceptual sub- 
tlety. In ordinary quantum mechanics, the steepest descent approximation provides the 
leading term in the transition amplitude in an h expansion. In our case, the action S that 
features in the path integral (3.15) itself depends on h through £o ~ y^^y^HG, while the 
h expansion of Appendix A assumes that the action does not change as h tends to zero. 
Therefore, to directly apply the result of Appendix A, now we have to take the hmit — > 
while keeping fixed. Hence we will obtain the leading term in the extraction amplitude 
in the approximation /I — )■ 0, 7 — )■ oo keeping 7^/1 fixed. To emphasize this subtlety, we 
will use inverted commas, as in 'classical limit' and 'semi-classical approximation' while re- 
ferring to this limit. Let us briefly explore the meaning of this limit. In classical general 
relativity, 7 — > 00 corresponds to ignoring the new term in the Hoist action for general 
relativity, in comparison with the standard Palatini term. What about the 'semi-classical' 
approximation? Eigenvalues of the volume operator are given by {^TiGloh)n where n is a 
non-negative integer. Therefore, in the 'semi-classical limit' the spacing between consecutive 
eigenvalues goes to zero and v effectively becomes continuous as one would expect. Finally, 
states that are relevant in this limit have large n, just as quantum states of a rigid rotor 
that are relevant in the semi-classical limit have large j. 
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A. The Hamilton-Jacobi function S'lo 

To calculate the 5*10 term, we need to solve the equations of motion obtained from the 
action (3.16), then evaluate the action along those trajectories, and finally express the result 
in terms of the given initial and final points. The (positive frequency) trajectories which 
solve the equations of motion can be written in terms of two integration constants, i/b and 

0s, as 

1/(0) = i/BCOsh(Vl27rG'((/)-0B)), (4.2) 
6(0) ^ 2sign(^B) ^^^-i^^-VmG^4>-M)_ (4.3) 

These solutions have several interesting features. 

(i) As seen from the cosh dependence of the volume, these trajectories represent bouncing 
universes, with 0^ and ub giving the scalar field and volume values at the bounce point. The 
minimum volume vb is related with the scalar field momentum p by {vbI = 2£op/Vl27rG. 
Note that if ub is positive (resp. negative), then i/(0) remains positive (resp. negative) for 
all 0. For concreteness we will focus on trajectories with positive ub- 

(ii) z/(0) can vanish only on the trajectory with Ub = i.e., z/(0) = for all 0. Thus if 
we begin with the initial state Ui ^ 0, 0^, there is no (real) 'classical' trajectory at all with 
lyf — for any value of 0/. 

(iii) Whereas in general relativity all trajectories begin at the big-bang — they all tend to 
z/ = as — >■ — oo — it is obvious from (4.2) that all our trajectories tend to i/ — >■ oo in this 
limit (except for the trajectory z/(0) = V0). 

(iv) Recall that in full LQC, states which are sharply peaked at a low curvature configuration 
for large values of remain sharply peaked on certain 'effective trajectories' for all [5]. 
These are among solutions (4.3). 

(v) The relation between u and given in Eq. (4.2) coincides with the expression for the 
expectation value of the volume operator at a given scalar field value in any quantum 
state of LQC [21]. 

Evaluation of the action along these solutions can be greatly simplified if one integrates 
by parts the term — dr ^bu in (3.16). Then, using the equations of motion, the terms ^bu 
and p0 cancel each other and the action evaluated along the solutions is just given by only 
the boundary term, 

S\o^l {l^^b^ - Ufbf) . (4.4) 

To express (Sjo in terms of initial and final configuration variables, we need to solve for the 
constants Ub and 0^ in terms of Uf, (f)f; Ui, (f)i. Without loss of generality we can take 0,j = 
and 4>f = ifi (by setting (/? = 0/ — 0j at the end, one recovers the general case). Then we are 
led to solve the equations 



Ui — vb cosh(— V 127rG0B) 
— i/BCOsh(-\/l27rG(</7 — 0b)), 



(4.5) 
(4.6) 
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for vb and 0b in terms of the initial and final data: 



(4.7) 



i/B = — (_e-Vi2^v' + ^). (4.8) 



sinh(Vl27rG'(^)| V V ^> 



Clearly, i'b,4>b are real for any given initial configuration {ui^cpi) if and only if the final 
configuration satisfies 

g-yi2^|<^| < ^ < gyi2^|.^|_ (4 9) 

This is the necessary and sufficient condition for the existence of real trajectories. Let us 
first focus on the 'classically' allowed region (4.9). Using (4.3) to express hi and 6/ appearing 
in (4.4) in terms of the initial and final data (i//, 0/; i^,, 0j) we obtain the desired expression 
of the Hamilton- Jacobi function: 




(4.10) 

where (p = (f)f — (pi. The 'classically' allowed region consists of the upper and lower quarters 
in Fig. 1. For Vf^ipf in these two quarters, 5*0 is real and thus the amplitude (4.1) has an 
oscillatory behavior. Outside these regions the action becomes imaginary and one gets an 
exponentially suppressed amplitude. Thus, the situation is analogous to that in quantum 
mechanics. 

For completeness, let us now discuss the case where the final point lies in the 'classically' 
forbidden region, this is to say the situation where the boundary data satisfy 

^ ^ ^-Vvi^M or ^ > e^^^l^l. (4.11) 

For concreteness let z/j be positive as in Fig.l but now there is no restriction on the sign 
oi Vf. To find extrema of the action that join the initial and final configurations satisfying 
(4.11), we can follow the semi-classical procedure used to calculate tunnehng amphtudes in 
familiar systems and allow paths with imaginary momenta. Let us define 

b = ib, p = ip, d = ia, S = iS. (4-12) 

Eq. (3.16) then implies 

S^Jdrl^p^- hu + &(^p'- ^ttGu^'^^^P^^ j . (4.13) 

We now consider the case when the tilde quantities are real and compute the stationary 
trajectories of S. The 'positive frequency' (i.e. p > 0) trajectories are parameterized by two 
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FIG. 1: For fixed {vi,4>i), the (dashed) curves vj = i/-e^^^^^^(<t>f-<t>i) divide the (z//,(/)/) plane into 
four regions. For a final point in the upper or lower quarter, there always exists a real trajectory 
joining the given initial and final points (as exemplified by the thick line). If the final point lies 
on the left or right quarter, there is no real solution matching the two points. The action becomes 
imaginary and one gets an exponentially suppressed amplitude. 



integrations constants, Uo and 0o, and take the form 



z/(0) = z/oSinh(Vl27rG(0-0„)), (4.14) 
5(0) = 2sign(z/,) tanh-^(e-^^l'^-'^°l). (4.15) 

They represent universes that go through a singularity at = (po, where the volume vanishes 
and b diverges. As in the 'classically' allowed region we have |z/o| = 2£op/VT27rG. In terms 
of the initial and final data, the integration constants are 



,v^.„ = h^-'^f/-^ (4.16) 



sign(t// - 
sinh(Vl2vrG(/7) 




(4.17) 



which, as expected, take real values in the 'forbidden' region (4.11). The action can then 
be evaluated as before. Although now the paths encounter a divergence in b, the integral 
(4.13) is convergent and given by the tilde version of (4.4). (Moreover, the product z/(0)6(0) 
is always finite and vanishes at = 0o.) For the case when ip > and z/j/z/j < e-v'i27rG</3 ^^^e 
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result is 

(4.18) 

Similar expressions hold for other regions. For instance, in the </? > 0, Vf/vi> e^'^^^^f case 
the action takes the same form, except that the arguments of the tanh~^ functions are the 
reciprocals of the ones appearing in (4.18). 

In all cases, ^lo is negative; the extraction amplitude is exponentially suppressed for 
paths in the 'classically forbidden' regions. But as we approach the dashed curves marking 
the boundary of the 'classically' allowed and forbidden regions, the action S goes to zero. 
In particular then, from Fig.l it may appear that for any given z/j there is a significant 
probability of reaching the singularity i/y^ = for large ip = (pf — (f)i. However, as is common 
in more familiar systems, the steepest descent approximation also becomes poor in a neigh- 
borhood of the dashed curves! Indeed, wc know from full LQC that (in the deparameterized 

framework) the expectation value of tends to infinity for large ip. More generally, plots 
of the exact extraction amplitudes in Sec. IV C will show that the amplitude is always sup- 
pressed in the classically forbidden regions. Thus, while the steepest descent approximation 
provides much physical insight, it is by no means a substitute for the full quantum theory. 



B. det^^^lo and the WKB approximation 

To compute the amplitude (det 5'^S\o) one would need to perform some regularization 
in order to deal with the infinite dimensional determinant. This is can be done in ordinary 
quantum mechanics or field theory, and should as well be doable here. We will however take 
a different route and calculate this factor by means of the WKB approximation [26, 27]. 

Note that the extraction amplitude A^Uf, (pf] z/j, (pi) can be thought of as a physical state 
if one takes the initial data as fixed parameters and the final data as arguments of the 
wavefunction: the family of states 

^•^.A'^f, <l>f) ■= ^{^h 0/; ^i. 00, (4-19) 
parameterized by Vi and 0j satisfy the constraint equation 

^*..A = 0. (4.20) 

The h expansions underlying the desired WKB approximation are discussed in Appendix A. 
We begin with the ansatz for the physical state: 

*..a(^/,<^/) = a(^/,0/;^.>.)e^^(^^'*^^'^-'^'^^''^'^- (4.21) 

Following the procedure of Appendix A, the imposition of the constraint equation (4.20) to 
zeroth and first order in h leads to the following equations for a and W: 



C{i'f,(pf,di,fW,d^fW) = 0, and Cx a = 



(4.22) 
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where 

C{iyf,(f)f,bf,pf) ^p} 
is the 'effective constraint', and 

A 

is the vector field on configuration space qf = {uf, (pf) obtained from the Hamiltonian vector 
field of the constraint. 

The first equation is the Hamilton- Jacobi equation and, as expected, one can check that 
^lo given by (4.10) solves it. The amplitude a is determined by the second equation together 
with the condition 

a{ui, (pi] Uf, (pf) = a{uf, (pf, Ui, (pf) (4.23) 
which follows from the fact that ^ui,<f>i{j^f,<Pf) = '^i'f,<i>f{^i,(pi)- 

a = \uf (e-^^ - ^) (-e-^^^^ + ^)|-V4. (4.24) 

This is the factor we identify with {det5^S\o)~^^^ . Note that this quantity diverges at 
Vf — i/jC^^'^^^^'^ (dashed lines in Fig. 1) where the amplitude goes from oscillatory to 
exponential decay behavior. Thus, the WKB approximation can be valid only away from 
the dashed lines. This simply mirrors what happens in the WKB approximation in ordinary 
quantum mechanics. 

To summarize, we have succeeded in finding a saddle point approximation of the path 
integral as in equation (4.1). The determinant factor was not calculated directly but by 
matching with the terms of a WKB expansion. Therefore, we will call the resulting approx- 
imate extraction amplitude ^4^^^: 

>iw™(i.^,0^;i.„0,):=aet^l°, (4.25) 

where a is given by Eq. (4.24), ^lo by Eq. (4.10) and as before (f = (pj — cpi. We now proceed 
to numerically compare this approximate amplitude with the exact one. 

C. Comparison with exact solution 

One of the advantages of the model under study is its solvability [21]. In particular, it 
is possible to obtain a closed form expression of the extraction amplitude A{i'f,4>f;i'i,4>i) 
[15]. This is displayed in Appendix B. We calculated the exact solution numerically and 
compared it with the saddle point approximation obtained in Sec. IV B. We found that there 
is a good agreement away from the dashed lines of Fig. 1 which mark the transition between 
the 'classically' allowed region to the 'classically' forbidden one. Along the dashed line, 
however, the WKB amplitude diverges and the approximation fails badly just as in ordinary 
quantum mechanics. 

We illustrate these results in figures 2, 3 and 4. In the first two figures we plot of the 
real parts of the exact and WKB amplitudes as a function of Uf and 0/ respectively, for 



- ^ttG — Uf 



X = 



dC 



dpf 



15 



;i,=10; yj ilTiG iPf= Log(5) 

Re(A) 




FIG. 2: Real parts of the exact and WKB amplitudes are plotted as a function of the final volume 
rif := Uf/Mofi- The exact amplitude (dots) has support on the 'lattice' {vf — i/j)/4£o^ S At 
rif = fi-Q^'^'^'^G(<l>f-(f>i) _ ^]2ere is the transition from oscillatory to exponential behavior, and the 
WKB amplitude (solid line) diverges. (It also diverges at Uf = rii e~^i27rG(0/-9ii) _ 2/) Here, (/>i,fi 
and (pf are kept fixed: (pi = 0, := Vi/Moh = 10 and ^\2nGcj)f = log 5. 

«,=10; ;v=20 

Re(A) 



-0.1 : 
-0.2 : 




FIG. 3: Real parts of the exact and WKB amplitudes are plotted as a function of the final scalar field 
(j)f (dots and solid lines, respectively). Here (pi^Vi and uj are kept fixed: (pi = 0,ni := Ui/Moh = 10 
and rif := Vf/Moh = 20. The WKB solution is the curve diverging at \/l2iTG(f>f = log(z^y/fj) = 
log 2. 



fixed values of the remaining variables. The exact amplitude shows a sudden transition 
from oscillatory to decaying behavior. If one had access only to the exact result, this 
behavior would have seemed rather puzzling. The WKB approximation provides a physical 
understanding of this behavior. Thus, not only does the WKB approximation reproduce 
the qualitative behavior of the exact extraction amplitude away from the dashed lines of 
Fig.l, but it anticipates that the dashed lines mark a boundary between two quite different 
behaviors of the exact answer and provides a physical understanding of this difference. 

What can we say regarding the regime of validity of the saddle point approximation? 
From the path integral perspective, we expect it to be valid whenever S/h^ 1. From Eq. 
(4.10), we see that 5*10 scales with the volume times a coefficient which can be interpreted 
as measuring the departure from the dashed lines Uf = i/je^^"^^'^*^''^'. At these lines 5*10 = 
and, as we have just seen, the approximation totally breaks down. As we depart from these 
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FIG. 4: Comparison of the exact and WKB amplitudes for initial and final configurations with 
Uf = 2ui as a function of rii := Ui/Mofi- As rij increases, the WKB solution (continuous line) 
becomes closer to the exact solution (dots). (The distance between the amplitudes oscillates, but 
overall it decreases.) This calculation was done for = and ^/ = l/y/12-nG. 



lines, 5*10 takes a nonzero value, and its scale is given by the initial and final volume. For 
instance, if we keep the ratio Vf/vi fixed, the action grows linearly with z/j and we expect 
the approximation to improve as z/j increases. This behavior is indeed observed, an example 
of which is display in Figure 4. Thus, the standard expectations on the validity of the WKB 
approximation are all borne out. 

V. DISCUSSION 

For constrained quantum systems, the extraction amplitude encodes full quantum dynam- 
ics: It enables one to extract physical quantum states from (suitably regular) kinematical 
ones and define the physical scalar product between them. In this paper we considered a 
solvable model in LQC. Following Feynman, we began with the expression of the extraction 
amplitude A{vf, 0/; i/j, 0j) in the Hilbert space framework and obtained two equivalent forms 
of a phase space path integral for it. In the first, one is led to integrate over paths that are 
very different from those one might have expected from the Wheeler-DeWitt theory: The 
integral is taken over quantum paths defined by the spectrum of the operators related to 
the phase space variables. However, by a clever trick from the quantum theory of a particle 
on a circle [23], this path integral could be reduced to one over standard classical paths. 
This form is better suited for semiclassical considerations and is also closer to that used 
in the Wheeler-DeWitt theory. However in neither form is the weight associated with a 
path given by the Einstein-Hilbert action. Instead, it is given by an action that includes 
quantum geometry corrections. This is the key difference from the Wheeler-DeWitt theory. 
Note that such a change of action in the transition from the Hamiltonian quantum theory 
to a path integral can occur already in much simpler systems. For example, for a particle 
moving on a Riemannian manifold, dynamics generated by the standard Hamiltonian oper- 
ator, H = -{h^/2m) g^^VaVb, is correctly captured in the path integral framework only if 
one adds to the classical action an h dependent term that depends on the scalar curvature 
of the Riemannian metric [25]. 
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Prom the path integral perspective, these differences are directly responsible for the ab- 
sence of singularities in LQC, i.e., for the resolution of the apparent tension we began our 
discussion with in Sec.I. In the first form of the path integral, the classically singular paths 
are not contained in the range of integration since the paths in b are bounded above while in 
the classical singularity b = oo. Furthermore the action is no longer the classical action. In 
the second form of the path integral, the domain of integration does include singular paths. 
But now the fact that the action is corrected by quantum geometry effects becomes crucial. 
Indeed, in this case the equations of motion can be obtained explicitly by varying the action. 
These equations and their solutions describe bouncing cosmologies which are characteristic 
of the singularity resolution in LQC. Thus the exact results on singularity resolution in 
LQC are in complete harmony with the path integral intuition, once one realizes that the 
action that descends from the Hamiltonian theory includes quantum geometry corrections. 
Furthermore, because we have an additional constant in the theory — the Barbero-Immirzi 
parameter — it is meaningful to consider h expansions while retaining quantum geometry 
effects. This is achieved by sharpening the precise manner in which the limit is taken: 
^ — )■ 0, 7 — )■ oo such that = const. This h expansion enables us to introduce the 
WKB approximation which helps one understand features of the exact amplitude, e.g., the 
oscillatory versus damping behavior that occurs as one varies the final configuration {i/f, (pf) 
keeping the initial configuration (z/j, (pi) fixed. It also provides an 'explanation' of the sur- 
prising effectiveness of the effective equations [28] in LQC from a path integral framework. 

Thus, from the LQC perspective, it would be incorrect to simply define the theory start- 
ing with smooth metrics and matter fields and assigning to each path the weight that comes 
from the Einstein Hilbert action because this procedure completely ignores the quantum 
nature of the underlying Riemannian geometry. For a satisfactory treatment of ultraviolet 
issues such as the singularity resolution, it is crucial that the calculation retains appropriate 
memory of this quantum nature. This viewpoint can be traced back to full LQG and spin- 
foam models. In full LQG, quantum geometry is an essential feature already of kinematics. 
It is then not surprising that in spinfoams the histories that one sums over are quantum 
geometries (captured in appropriate, colored 2-compIcxcs). The weight that is assigned to 
each history is motivated by the Einstein-Hilbert action but does not descend directly from 
it, e.g. via a discretization procedure. Rather, one begins with a the action of a constrained 
BF theory which is classically equivalent with the Einstein Hilbert action but then incor- 
porates the (simplicity) constraint using considerations from quantum geometry (quantum 
tetrahedron, representation theory and interpretation of the Casimirs as eigenvalues of ge- 
ometric operators). Thus, the situation is parallel to the first form of the path integral we 
obtained in this paper. In LQC we were fortunate in that the integral over quantum paths 
could be recast into an integral over all paths in the classical phase space. This enabled us 
to carry out the steepest descent approximation and develop physical intuition for the qual- 
itative properties of the exact extraction amplitude. A similar reformulation of spinfoams of 
the full theory appears to be difficult. But if it could somehow be achieved, one would have 
a powerful tool both to probe semi-classical aspects of full quantum gravity and to develop 
valuable intuition for the ultraviolet properties of the theory. In particular, the resulting 
quantum geometry corrections to the full Einstein-Hilbert action would bring the difference 
between spin foams and perturbative path integrals into sharp focus. 
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Appendix A: WKB approximation for constrained systems 

Let us consider a system with phase space M^" and a single constraint C{q,p) = 0, to be 
thought of as the Hamiltonian constraint. We will assume that the C{q,p) can be written as 
a Taylor expansion in the p, as is the case for a large class of physically interesting systems. 

The kinematic Hilbcrt space is L^(]R") of normalizable wave functions \E'(g) = 
'i/{q^, . . . ,q"'). The elementary operators are the usual q^ and pj which act by multipli- 
cation and derivation respectively: 

q^xif{q) = gJ^(g) (Al) 

pMq) = -^^|^(?)- (A2) 

The physical states are then solutions to the quantum constraint: 

C* = 0, (A3) 

where C is an operator analog of C{q,p), obtained by replacing q,p with q^ and pj with 
a suitable choice of factor ordering. As is standard in the group averaging procedure, we 
will assume that C is self-adjoint. For unconstrained systems, the WKB ansatz provides 
approximate solutions to the Schrodinger equation. In this Appendix we will extend that 
method to obtain approximate solutions of (A3) where, again, h plays the role of small 
parameter governing the expansion. As one might expect, the main idea is to write both C 
and ^ in (A3) as expansions in h, and to collect terms having the same h power. Let us 
now explicitly calculate the zeroth and first order terms. 

The construction is as follows. First, the constraint operator C is written as a sum of 
'normal ordered' operators, in which all g's appear to the left of the p's: 

oo 

5 = ^(;iA)"c„(y ). (A4) 

n=0 

Here the C^s are functions on the classical phase space — for instance, Cq will typically 
be the classical constraint function — which are now 'evaluated' on the operators q and p 
according to the 'normal ordered' prescription indicated by the superscripts. Second, the 
unknown state ^(?) is written as the exponential 

^(g) = et^(«) (A5) 
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where the exponent is written as a power series in h: 

oo 

'5(?) = E(W'5n(g). 



(A6) 



n=0 



Since C{q,p) is assumed to admit a Taylor expansion in the p, so do Cn{q,p)- Imposition 
of the quantum constraint (A3) now leads to the following zeroth and first order equations: 



1 d'Co 



2 dpidpj 



P=dqSo 



dq'dq^ dpi 



P=dqS{) 



dq^ 



+ Ci(?,9,>So) = 



(A7) 
(A8) 



The zeroth order equation (A7) can be recognized as the Hamilton- Jacobi equation. The 
first order one, (A8), can be rewritten as follows. If we use the fact that C is self-adjoint, 
the condition = C, when applied to (A4) implies 



1 a^Cn 



2 dqWpj ' 



(A9) 



Using (A9), and writing a{q) :— e'^^^^^ Eq. (A8) can be written as a derivative of the 
function a along the vector field 



X 



dpj 



p=dqSo 



_d_ 

dqi 



as 



X{a) + -adivX = 0. 
2 



(AlO) 



(All) 



The divergence term in (All) suggest one to interpret a and ^ ~ aeft'^° as half densities 
on R". Then (All) is just the Lie derivative of a along X: 



Eq. (AS) <^ 

These are the equations used in Sec. IV B. 



Cxa = 0. 



(A12) 



Appendix B: Exact Amplitude 

In Sec.IV C we compared the exact extraction amplitude with the WKB approximation. 
In this Appendix we recall from [15] the expression that was used in the numerical evaluation 
of the exact amplitude. 

The first step in the calculation is to find the eigenvectors of G. They are given by |A;±), 
with A; > 0, satisfying the eigenvalue equation 



e|A;±) = l2T:Gk'^h^\k±). 



(Bl) 
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These vectors are not normalized; the decomposition of the identity reads 

I = /o°^2^™T(^ {\k+){k+\ + \k-){k-\). (B2) 

In terms of the 'volume basis' \Anioh) used in the main body of this paper, the vectors \k±) 
are given by 

{An£^h\k±) = I VW¥^kP±nik) ±n ^ ^ ^23^ 
where P„(A;) is the following (2n — l)-degree polynomial in k: 

ik 2n ^ 2n-l 



1 f/^" 

Pn{k) 



ik(2n)\ ds'^'' 



(t^I ^Y\{tk + m-l). (B4) 



We are now ready to present the expression of the extraction amplitude. For this, it 
is convenient to work in the deparameterized framework. In [15] it was shown that the 
extraction amplitude in the timeless framework coincides with the transition amplitude of 
the deparameterized theory: 

A{uf,(l)f,u,,(/)i) = {uf\e^^^^f'f-^'^\ui). (B5) 

Let us take z/j = Mohrii and i/f — Alghnf with rij and Uf positive integers, and define 
t := v^T2^(0/-0i). 

By inserting the complete basis (B2) in the right hand side of (B5), we obtain 

A(uf, 4^f- V,, = -2V^/o~ii4^ Pnsik) PnXk) ke''' (B6) 
= -2,/Tm]P^^i-idt)PnX-idt)J7^^)ke'''' (B7) 



,/n-nJ 



PnA-idt) PnA-idt) ¥'\l/2 - t—) (B8) 



where i/j^^^z) = dlogT{z)/dz, and T{z) = [z — 1)! is the Gamma function. 

This last expression (B8) was the one used to numerically compute the exact extraction 
amplitude for the plots in Sec. IV C. 
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